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We explore extensions of hybrid inflationary models in the context of supersymmetric D-term in- 
flation. We point out that a large variety of inflationary scenarios can be encountered when the fleld 
content is extended. It is not only possible to get curvaton type models but also scenarios in which 
different flelds, with nontrivial statistical properties, contribute to the primordial curvature fluctua- 
tions. We explore more particularly the parameter space of these multiple fleld inflationary models. 
It is shown that there exists a large domain in which signiflcant primordial non-Gaussianities can 
be produced while preserving a scale free power spectrum for the metric fluctuations. In particular 
we explicitly compute the expected bi- and trispectrum for such models and compared the results 
to the current and expected observational constraints. It is shown that it is necessary to use both 
the bi- and tri-spectra of CMB anisotropies to efficiently reduce their parameter space. 
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I. INTRODUCTION 

The new generations of high precision Cosmic Mi- 
crowave Background (CMB) measurements opens a new 
window to the physics of the early universe and more 
specifically to the physics of inflation. Up to now, the 
shape of the power spectrum has been the prime object of 
interest in order to confront theory and observation. Any 
observed departure from a scale free, Harrison-Zel'dovich 
type, power spectrum would indeed give a precious hint 
on the shape of the inflaton potential and subsequently 
on models of inflation [36]. But needless is to say that 
detection of primordial non-Gaussianities in the CMB 
temperature anisotropies or polarization would also carry 
precious indications regarding the inflaton sector. 

It has been shown however that in its simplest for- 
mulation inflation does not lead to significant non- 
gaussianities (NG). This result has been rigorously es- 
tablished at the level of the bispectrum over the last 
few years [HJl, J3| and recently it has been extended to 
trispectrum [j, |5|- To be more precise the low level of 
primordial NG that standard infiation is expected to in- 
duce is unlikely to be detected in practice due to the 
subsequent evolution of the metric whether it is at su- 
perhorizon scale, during or after the inflationary phase 
or during the complex stages of the recombination 0] ■ 

It appears then useful to have at our disposal sets of 
models in which significant primordial non-Gaussianities 
can be produced, that is an amount of NG that clearly 
exceeds what the nonlinear gravitational evolution is in- 
ducing. This appears to be the case for instance for mod- 
els built out of actions close to the Dirac-Born-Infeld one 
Q . Those models however rely on non trivial exten- 
sion of the kinetic part of the Lagrangian. Although that 
would be a ground breaking discovery to see such effects, 
it remains that the appearance of primordial NG does 
not require the use of such non-standard extension of 
the low-energy field equation. The curvaton model of- 
fers another path to possible large amount of primordial 
NG p^ . Here primordial non-Gaussianities are induced 



lately, at the time the energy density of the curvaton field 
overcome the inflaton decay products. The efficiency of 
such a mechanism therefore depends on sectors of the 
theory that are not directly relevant to the inflation. 

A third possibility has recently been advocated which 
is entirely determined in the inflationary sector. It is pro- 
vided by some models of multiple field inflations. The 
mechanism at play has been described in various papers 
pH [T^ . The difficulty so far is to identify more precisely 
models in which such a scenario could take place. In 
explicit Lagrangian (with canonical kinetic terms) were 
proposed for which significant non-Gaussianities can be 
generated during the inflationary phase. Although it was 
demonstrated that it was possible to get them during 
the whole inflationary stage, the most efficient (and in 
some sense most natural) model in which such a possi- 
bility was achieved was based on an hybrid inflation type 
construction. In this case significant non-Gaussianities 
can be transferred to the metric fluctuations at the very 
end of inflation, at the time of the tachyonic instabil- 
ity ^37j . The model proposed in [l^ was however purely 
phenomenological and was not motivated by high energy 
physics constructions. 

We propose here to explore a bit further high-energy 
physics motivations for multiple-field hybrid inflation. 
Hybrid inflationary models introduced by A. Linde fT3| 
provided a way to circumvent one fundamental issue in 
standard inflation with simple potentials, that is that the 
Vacuum Expectation Values (VEVs) of the inflaton field 
is bound to be of the order of the Planck mass during 
the inflationary period p^ . This is the case in partic- 
ular with quadratic or quartic potentials. With hybrid 
inflation the end of inflation is triggered by a tachyonic 
instability and the VEVs of the fields can be made as 
small as one wishes (at the expense of a power spectrum 
index very close to unity). Hybrid inflation has received 
some further theoretical justification because it naturally 
appears in various high energy setups. In the context of 
global supersymmetry, it is indeed possible to construct 
potentials which allow a viable inflationary period, either 
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based on the F- or the D-term of the potentials [l^, |T3| ■ 
The aim of this paper is to take advantage of such theo- 
retical frameworks to build models of multiple field infla- 
tion and explore their phenomenological consequences. 

The plan of the paper is the following. We first suc- 
cinctly review the construction and constraints associ- 
ated to supersymmetric hybrid models. We then present 
possible extensions of these models and their phenomeno- 
logical consequences. The amount of non-Gaussianity 
they can lead to is precisely computed and confronted to 
observational constraints. 



II. SUPERSYMMETRIC HYBRID INFLATION 

In a supersymmetric context, two scenarios have been 
identified that lead to infiationary models that reproduce 
hybrid type potentials. Those scenarios are the F- and 
D- term inflationary models [Tsj. They allow in partic- 
ular to have inflation while the field VEVs remain small 
compared to the Planck mass. In such a case the expres- 
sion of the low energy Lagrangian is safe from corrective 
terms originating from supergravity effects (e.g. through 
the Kahler potential) or high-order curvature terms. This 
is the domain we will stick to throughout this paper. 

To be more specific, the F- and D-term infiation are 
based on the specific global supersymmetry breaking 
mechanisms. Once symmetry breaking occurs, the po- 
tential has to be corrected for one-loop radiative correc- 
tion terms. These corrections are encapsulated in the 
the formula of Coleman- Weinberg derived in [l^. For 
soft supersymmetry breaking schemes (where the trace 
of the mass matrices remain unchanged) one is left with 
logarithmic corrections to the potential. This is precisely 
what we have for F- and D-term inflationary models. 

The F-term inflation is based on a superpotential of 
the form UM. 



W = \S(t>(t> - fi^S 



(1) 



where S, (j) and (j) correspond to the scalar degrees of free- 
dom. The resulting potential is the following. It comes 
from F-tcrm contribution only, 



V = X^\S\ 



\X(f)(f> - ^1- 



2|2 



(2) 



The energetic landscape it leads to shows that there ex- 
ists a flat direction along which the potential is constant 
(at tree order) but non-vanishing. It is for a such field 
configuration that the inflationary phase can take place. 
This flat direction is obtained for (p ~ (j> — 0. It corre- 
sponds to a stable minimum (e.g. all masses are positive) 
as long as 



(3) 



that is |5| > Sc with Sc = n/VX. Such a vacuum state 
obviously breaks global supersymmetry. It implies that 



corrective terms due to non-vanishing radiative correc- 
tions have to be added to the potential. At one-loop 
order, they approximatively read, 



loop 



^logM 
167r2 Sc 



(4) 



This is thanks to this extra contribution that such a set- 
ting can lead to a viable inflationary model. It induces 
a (slow) roll of the field S towards the origin. In this 
scenario inflation terminates [s^ when |iS| = Sc due to 
the tachyonic instability that then appears. During the 
inflationary phase the potential can be approximated by. 



VF-term(5) = X'S^ 1 



A2 



■log 



151 



(5) 



to a very good approximation. 

The mechanism that leads to a D-term inflation is very 
similar [17]. This model is based on the introduction 
of a U{1) symmetry group and a non- vanishing Fayet- 
Iliopoulos term leading to a non-zero D-term in the po- 
tential. Another family of inflationary models can then 
be constructed with the help of three flelds, one, S, with 
zero charge, a fleld (f) with a positive charge and a field cf) 
with a negative one under the U{1) symmetry and with 
the following superpotentiel, 

W = A5#. (6) 

It leads to the potential, 

V = Vf + Vd 

= x'\s\' {\w + \<f>n + >^'\M' + Y - 

In a way similar to the F-term inflation, the potential 
exhibits a flat stable direction when, 



- 0, = 0, \S\ > Sc (8) 



with 



5c = fVC. (9) 
In this case the inflaton potential is given by. 



v,_..„,,5,.:^(i + |l,„g|lj. (10) 

This is the context in which we will work in the following. 

A. Conditions to have a valid D-term inflation 

A number of conditions should be met for such models 
to provide us with valid inflationary models: the number 
of efolds it leads to should be large enough and the am- 
plitude of the metric fluctuations is constrained. These 
are obviously constraints that are generic for all infla- 
tionary models. There are however two difficulties that 
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one wants to avoid in the context of D-term inflation. 
Namely one does not want the VEVs of the fields to ex- 
ceed the Planck scale so that global supersymmetry re- 
mains a valid description of the Lagrangian'so'l and the 
formation of too massive cosmic strings at the end of the 
inflationary phase should also be avoided. 

The first requirement leads us to assume that not only 
Sc is small but also that the first term of the inflaton 
potential, the constant one, is dominating over the second 
contribution. This is possible if 



5' « 1- 



(11) 



With such a constraint the end of inflation takes place 
when |iS| = Sc- The number of efolds between horizon 
crossing for the modes of interest (the ones that are re- 
sponsible for the observed large-scale structure) and the 
end of inflation can then be easily computed. It is given 
by 



N, 



87r2 



(12) 



PI. 



where 5* is the modulus of S during horizon crossings 
and Mpi is the Planck mass, M|,j = l/(87rG'). The 
amplitude of the metric fluctuations TZ* for the modes of 
interest is observed to be of the order of 10"'' (see next 
section for more details). It implies that 



47r2 £_S, 



10- 



\/3 ffAff,! 

These two relations lead to the constraint, 



3 M4j^ 



87r2 



^'^^Pl. 



10' 



(13) 



(14) 



C is fixed, of the order of IQ-^ M|,j 

e/(A^A%) 



This expression exhibits two regimes. One in which 
^/{X^M"^^ ) < iVe/(87r2) and then the energy scale of 

and one in which 
> and where ^ ~ 0.5 A^/^ IQ-^ Af|,^ . 
But ^ is precisely the energy scale of the strings that form 
at the end of infiation. The constraint on the string con- 
tent of the Universe favors the second case (see [20, HH 
for details. It is to be noted also that it might be possible 
to circumvent this constraint). 

It is to be noted that in this latter case 5* is nearly 
equal to Sc- More precisely we have. 



5c 1 



167r2 52 



(15) 



We will see that this is of importance for the phenomeno- 
logical consequences of D-term extended models. 



III. EXTENSION OF THE FIELD CONTENT 

Here we now propose to extend the field content in- 
volved in the inflationary sector of the theory. Following 



a conservative approach regarding the high energy sec- 
tor of the theory wc allow ourselves to introduce up to 
cubic terms only in the superpotential. For instance a 
superpotenticl of the form. 



W = Ai5i^ 



(16) 



is a priori legitimate. It involves two light fields that 
coexist during the inflationary phase. If ^2 is smaller 
than the Hubble value during the inflationary phase then 
super-Hubble fluctuations can be generated in the 52 
field. Whether those fluctuations will be observable or 
not depends on the subsequent evolution of the 52 fleld. 
In particular, because ^2 is a massive fleld, its energy 
density can eventually dominate over the inflaton de- 
cay products if those are relativistic particles. The ^2 
field can then imprint its fluctuations on the metric fluc- 
tuations before it finally decays. This is the curvaton 
mechanism. This possibility, and its phenomenological 
consequences, have been extensively described in the lit- 
erature [13] ■ They depend largely on the couphng of the 
extra fields to the other fields of the theory. Such a theory 
is then not entirely predictive from the sole knowledge of 
the inflationary sector of the theory. 



A. Multiple field models 

The expression is however not the only possible 
extension of ([5]). For instance nothing prevents the intro- 
duction of multiple light fields that are coupled together 
to the same charged U{1) fields, 



W ■ 



(17) 



Obviously if Vi is small enough, the corresponding Si field 
can participate in the inflaton (depending on its initial 
VEV) . The corresponding upper bound for Vi for such a 
possibility to occur is defined in such a way that, when, 
say, the VEV of Si is below the Planck scale, the con- 
tribution of the quartic potential it induces is negligible 
against the radiative correction term. It leads to the con- 
straint. 



(18) 



The fields for which Vi is above this bound will rapidly roll 
towards the origin but they still can develop significant 
super-Hubble fluctuations as long as Ui is smaller than 
unity. This would not be the case otherwise. 

It is then possible to distinguish three sets of fields, 
those with a very small coupling constant v\ they can 
potentially be part of the inflaton field; the ones with 
intermediate values ; they will not contribute to the in- 
flaton but still develop significant super-Hubble fluctu- 
ations and finally those with a large coupling constant 
will not develop any significant fluctuations and will not 
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play any role. In the following we will be interested in is one field Si with a vanishing coupling constant Pi and 

the second set of fields. one field ^2 with a large - but still smaller than unity - 

To simplify the presentation let us assume that there 1/2- 

I 



Then the potential takes the form, 

V = V^i-loop + ^^ |cos6i5i + sin 61521^ (101^ + |0|^) + X"^ cos"^ 9 \<j)f\^f + |j^25| + Asin^. 



^2 



(19) 

where 9 is the mixing angle of Si and 52 encoded in the at jjarameters 



As for the previous case, </> is a very massive field whose 
VEV is driven to 0. The previous expression then sim- 
plifies and we are left with the effective potential, 

V = Vi_iQQp + iyi\S2\^ + X^\cos9Si+am9S2f\W 

involving the (complex scalar) fields Si, S2 and (j>. We 
see that it involves a new term, the self coupling term 
of the 52 field. It drives the VEV of this field to dur- 
ing the inflationary period. The expression of Vj^-Joop 
then left unchanged compared to the single inflationary 
field case. As for the curvaton model, the 52 direction 
behaves like an isocurvature direction during the infla- 
tionary period. Because of the mixing term in the Si-cfi 
coupling though, the ^2 fluctuations can lead to visible 



effects irrespectively of the subsequent evolution of this 
field. To be more precise, the inflation terminates when 
cos 9Si + sin 9S2 reaches the critical value Sc- In the con- 
text we are considering, the inflationary period then ends 
almost instantaneously because of the tachyonic instabil- 
ity that subsequently develops [H, HI] . The end time of 
inflation is then modulated by both fluctuations in the 
Si and ^2 directions. This is the mechanism with which 
initial isocurvature fluctuations are transferred into the 
adiabatic modes. 

The metric fluctuations that are then induced are a 
combination of Si and ^2 fluctuations. As long as the 
only leading order metric effects are taken into account, 
the formal expression of the metric fluctuations can be 
easily computed using the 5N formalism [2^, [2^ . The 
induced metric fluctuations then read 



SN{5Si,SS2) = I . _ H{S[°'^ +SSi,6S2)dt- f ^ H{s[°\o)dt, 
^inflation traj. ^inflation trai. 



inflation traj. 



(21) 



where S['^'^ is the zero mode trajectory of the field. In the context we are interested in the expression of SN can easily 
be explicited at leading order in the field fluctuations. It is given by. 



SN 



Horizon crossing 



Sip - 



end of inflation 



taii9 Sxi 



(22) 



where ip and x are the (canonically defined) fluctuations 
of the fields Si and ^2 in the direction (in the complex 
plane) of 5^°^ . This direction, without loss of generality, 
can be defined as the real axis, 

ip^V2 Rc(5i), (23) 

and the real part of ^2, 

XI = \/2 Rc(52), and X2 = M^s). (24) 

The imaginary parts of those fields are genuine degree 
of freedom that will develop super-Hubble correlations 



as well. They will affect though the expression of SN at 
quadratic order only, as any other couplings to the metric 
would. At the level of our description the imaginary part 
of if will not play any role. The one of x will do however, 
because it affects the non-Gaussian properties of xi as 
we will discover. 

What is important here to realize is that in 
the parameter domain favored by the constraints 
on the cosmic strings contribution to the CMB 
anisotropics as discussed previously, the two coefficients, 

I horizon crossing ^^"^ '^-^^/^¥'lend of inflation 
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are, almost, equal. Moreover the two fields are indepen- 
dent of one another; their fluctuations have roughly the 
same spectrum, to the difference in their mass, e.g. 
parameter. The generated metric fluctuations are then 
the superposition of these two contributions the relative 
weight is driven by a free parameter, the mixing angle 
between the two fields, 9. The amplitude of the induced 
metric fluctuations per unit log scale in fc, is then 
typically given by [40|, 



These parameters however vanish as soon as g is small 
and we are left with a scale free Harrison-Zel'dovich spec- 
trum in all cases. 



B. The induced non-Gaussianities 



Vo-- {l + tan^ 



,1/2 3i/3 



2 cos 9 



(25) 



This is this number that is constrained by the observa- 
tions: Vq is of the order of 2 x 10~^ (see for instance 

One can be a bit more precise. The first term is ex- 
pected to be a random field of power spectrum of index. 



rti = 1 — 6e + 277, 
and the second of power spectrum index, 
n2 = l-2e + 277i, 



(26) 



(27) 



where e is the same in the two cases (it is due to the 
variation of i7). 



H 



(28) 



and rj are the masses of the fields in units of H in re- 
spectively the adiabatic direction and the transverse di- 
rection, 



V 

^'^Pl. 



V 



(29) 
(30) 



The resulting index, in first order in slow-roll parame- 
ter is then. 



1 + cos^ 9{2f]i ~ 6e) + sin^ 9{2f^2 - 2e) 



(31) 



The most interesting consequences of such a family of 
models comes from the fact that they can potentially 
induce significant non-Gaussianities. By significant we 
mean that such non-Gaussianities can be much larger 
than what generic single field inflation predicts. The 
amount of non-Gaussianities which is generated is ac- 
tually a direct consequence of the formula Indeed, 
the self couplings of the field x can be large. We are 
typically in a situation described in [l^ . Those intrinsic 
non-Gaussianities are not due to any coupling of the field 
to the metric but to the self-coupling this scalar degree 
of freedom in a quasi-de Sitter background. 

In particular, it is clear that 882 can develop signifi- 
cant high order correlation function, its four point can 
be computed explicitly. Finite volume effects though can 
induce as well non-vanishing three-point function. As 
shown in [T^ it amounts to shift the minimum of the 
field to a non zero value. The effective potential in the 
transverse directions is actually expected to be 



nXi,X2) = ^ [{Xi+X,f + {X2+X2?]\ (32) 



where Xi and X2 are in essence random variables but that 
can be considered as fixed for our observable universe. A 
classical stochastic approach can be applied to infer the 
probability distribution function of its value. The Fokker 
equation for the joint distribution of Xi and X2 can be 
derived from the joint evolution of Xi and X2- It is given 
by, 



dV_ 



1 

3H 



_d_ ( dV{xi,X2) ^ 
dXi \ dxi 



dX2 \ 9x2 



V 



(33) 



r 



where t is the physical time. The late time solution of for the correctly normalized distribution function. For 
this equation is the bispcctrum, wc will see that the parameter of interest 

is Xi- Its distribution function, once marginalized over 



l^(Xi,X2) = TTTT-T^exp 



, (34) 
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FIG. 1: Shape of the one-point distribution function of Xj for 
H = 1 and = 1. The solid line corresponds to the equation 
(|35[) where the imaginary part of x has been integrated out. 
The dashed line corresponds to the case where x is real. 



X2 is given by, 

Mxi\ 



nxi) 



Ki 



9 9 -4 



cxp 



? 2 -4 

xi 



(35) 

where Ki is the modified Bessel function of the second 

4 

kind of index 1/4. The shape of this PDF is shown on 
Fig. [T]for = 1 and = 1. 

For generic initial conditions, and for a large enough 
number of efolds prior to the horizon crossing of the 
modes of interests, Xi is expected to be distributed ac- 
cording to ([35]) . The excursion domain for Xi is therefore 
typically [—0.6H/^/l',0.6H/^/U\. As suggested before Xi 
can be seen as a free variable in a situation similar to that 
encountered in the curvaton scenario. Note that, from 
the potential ((32)) . the field xi acquires a mass equal to 
5i/2jy^^ It is generically small compared to H ensuring 
that the field can indeed develop super-Hubble fluctua- 
tions. 



C. The perturbative regime 

We are interested here in the high order (e.g. three- 
and four-point) correlation functions of the Fourier 
modes of the field x and consequently 5N. For a field 
cr(x) we define its Fourier transform CT(k) by, 



a(k) 



{2n) 



rcr(x) e 



— iks 



(36) 



The power spectrum P^lk) and the higher order spec- 
trum Pn{k) are defined by 



(cr(ki)cr(k2)) = (5Dirac(ki 
(CT(ki) . . . Cr(k„))c = (5Dirac(kl 



-k2)P'^(ki), (37) 
••• + k„)i^„"(ki,...,^) 



where c denotes the connected parts. 

A full quantum calculation can be done at tree order in 
the perturbative regime (see last section for details) . In a 
de Sitter background the exact shape can be obtained for 
the three and connected part of the four-point function 
for the potential ([5^ . In the super-Hubble regime they 
read gj, 



P3^^(ki,k2,k3) = - 



i/2 



[P'^i(fci)F'^i(fc2) +sym.; 



(39) 

where Ng is the number of efolds between horizon cross- 
ing and the end of inflation. 

As stressed in previous papers, what drives the ampli- 
tude of fluctuations is the value of v'^N^ if S is of the 
order of H . This result can be expressed phenomenolog- 
ically through an /nl parameter. It is defined in such 
a way that the bispectrum i-'3^(ki, k2, ka) of the metric 
fluctuation reads, 



P3^^^(ki,k2,k3) = 2/nl [P(fci)P(fc2) +sym.] , (40) 

when written in terms of the metric power spectrum. 
Relation p9|) implies that. 



/nl 



H 2Vo 



(41) 



From the adopted definition of /nl: the metric fluctu- 
ation will be significantly non-Gaussian if Vq /nl ap- 
proaches unity whereas standard inflationary physics im- 
plies that /nl is of the order of unity '3'| . 

Within the framework of this model, it is also possible 
to compute any higher order correlation functions. When 
Xi and X2 both taken into account, there are two 
terms contributing (and from a field theory point of view, 
it means that there are three- as well as four-leg vertices 
that contribute). 



They read. 



Pf (ki 



9x? 



xl 



9H^ 



P(fci)P(|ki+k2|)P(fc3)+sym.] 



2iy^ N,. 



This leads to 



pr(ki, 



,k4)=4/j 



^Xl 



\^xh 



[P(fci)P(|ki + k2|)P(fc3) + sym.] + GffNL [P(fci)P(fc2)P(fc3) + sym.; 



■ sym.J . 
(42) 

(43) 



k k k k K K 

pxi(ki,k2,k3,k4) = > — _> — <r^^y"^^x_ 

ki-^ .ka k ,k3 ''^•.-^ 

P^l(ki,k2-k3,k4) = J \* J \ ^ sym. + + sym. 

k* \k* *k4 k,* k> 

FIG. 2: Diagrammatic representation of the contributions to the connected four point functions of the xi component of the 
complex field x- The upper line show the Feynman type diagrams that have to be considered in a quantum field approach 
(for details see for instance [1^). For tree order calculations one can equivalently use a classical approach (provided the initial 
stochastic fields have the properties derived from the free field quantum solutions) . The bottom line shows the resulting diagrams 
where each vertex point then represents a given order in the initial field. The mode dependence of those terms can then easily 
be derived in the super- horizon limit: the three diagrams that are represented give respectively P^^ (k2)P^^ (l^i + \i.2\)P^^ (^4), 

P''l(fc2)P'^=(|ki +k2|)P^Kfc4), P^^{k2)P^^{k'i)P^^{ki). 



The two terms of this equation correspond to two differ- 
ent geometries. The coefheient ^nl reads, 



5NL = -V Ne 



sm 



(44) 



Note that the expression of the bispectrum and the sec- 
ond term of the tri-spectrum is what a nonUncar trans- 
form of the type, 



SN^5N + /nl {SNf + gNL {SNf 



(45) 



would have given. Such a transform would however lead 
to for the coefficient of the first term of the tri- 

spectrum. The difference is due to the fact that the re- 
sulting metric fluctuations are actually the superposition 
of different fields, not only the nonlinear transform of a 
single one. 

The relative weight of these two contributions depends 
in particular on the value of S, and therefore on the pe- 
culiar realization of our own universe. These results are 
valid only in the perturbative regime, e.g. with respect to 
the coupling constant. It appears that the expansion pa- 
rameter is different during sub-Hubble physics and dur- 
ing super-Hubble evolution. In the former case, z/^ is the 
parameter to be small. When it is large the field fluctua- 
tions simply vanish away. At super-Hubble scale however 
one finds that the effective expansion parameter is actu- 
ally v'^Ng which is obviously much larger than i/^. It 
limits the formal validity of the previous expressions. It 
also opens a new specific phenomenological domain the 
tentative description of which is the subject of the next 
section. 

In Fig. [21 we compare the computed amplitude of the 
bi- and tri-spectrum to the current (with WMAP) and 
expected (with Planck) constraints. The left panel cor- 
responds to a mixing angle 9 — 7r/4, that corresponds 
to an equal contribution of the two fields to the metric 
fluctuations, and the right panel to a small mixing angle. 



9 ~ 0.1. We restricted the parameter space to two vari- 
ables, v and Xii assuming X2 — ^ (it anyway contributed 
only weakly to the amplitude of the tri-spectrum) . The 
dotted lines show the location of the expected generic 
value for Xi- Values that would differ too much (in log- 
arithmic space) from this location would demand some 
fine tuning in the peculiar random realization of the uni- 
verse we live in. 

In both cases the thick solid line shows the current 
constraint provided by the amplitude of the bispectrum. 
It naturally limits a combination of Xi ^-nd v. When the 
mixing angle is large, left panel, it generically leads to 
small values of v. The use of a perturbation theory is 
then fully justified. When the mixing angle is small the 
observations are much mess efficient in constraining v. It 
can be as large as unity. It is then necessary to reconsider 
the calculations that have been done to take into account 
non-perturbative aspects of the super-Hubble evolution 
of the fields. 

The Planck mission is expected to provide us with 
much more stringent constraints not only on the bi- 
spectrum but also on the amplitude of the tri-spectrum. 
In Fig. [3]the dashed lines show the location of the bispec- 
trum constraint provided by Planck. Regarding the tri- 
spectrum we have not attempted to take into account the 
different geometrical patterns that appear in its theoret- 
ical expression. As a result the amplitude of the trispec- 
trum results of the summation of two terms of opposite 
signs. There is therefore a location in the parameter 
space where the tri-spectrum effectively vanishes. It is 
shown as a long dashed line. The constraint that obser- 
vations would provide is given in terms of tnl which is 
set to be equal to Gg^L + 4/nl/ sin^(6'). The gray area is 
the reg ion that Planck could be able to exclude according 
to [23] ■ It appears clearly that for low values of Xi , the 
tri-spectrum is more effectient in constraining the ampli- 
tude of the coupling constant v. In the rare event tails 
for Xii ^rid along the cancelation line, the bispectrum 
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FIG. 3: Exclusion diagrams for parameters v and x for 9 = n /4: (left panel) and for 9 = Q.l (right panel). The locations of 
the dotted lines where x is equal to its expected one a fluctuation. The gray areas and solid or short dashed lines correspond 
to the exclusion zones, obtained by WMAP (solid line) or expected by Planck (short dashed for bispectrum, gray areas for 
tri-spectrum). The bispectrum constraint corresponds to a straight line (of slope —2); the trispectrum is more complicated due 
to two competing terms in the trispectrum. The long dashed is the location where the terms cancel. We adopted the results 
of [2^ on the upperbounds the Planck mission is expected to provide, /nl = 5 and tnl ~ 560. 
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FIG. 4: Example of shapes of the one-point PDF of the local 
value of X in case it underwent a non-perturbative evolution 
at super-horizon scales. The plots correspond to parameter 
values vNe ~ 1.5; the long dashed line to x = 0; the solid 
hne to Xi = 0.5, X2 ~ the short dashed line to Xi = 0, 
X2 = 2 (values of x are given in units of H). The thin solid 
line is a Gaussian distribution of similar width. The resulting 
PDF would be the convolution of one of the first distribution 
with a Gaussian one with arbitrary relative amplitude. 



is more efficient. The two observations appear therefore 
very complementary. For low values of 9 though, it is 
obviously more difficult to get strong constraints of the 
coupling constant u. In this case the field S2 could ac- 
tually induce primordial NG when it enters a classical 
non-perturbative regime. 



D. The non-perturbative regime 

The aim of this section is to point to the existence of a 
regime where the statistical properties of the field ^2 are 
transformed due to its super- Hubble nonlinear evolution. 
Such properties could be transferred in the observed met- 



ric perturbation in case of a small mixing angle. 

Let us recall that at super-Hubble scales the evolution 
equation for each component of the field is given by, 



(46) 



for each component of the field and as long as the slow 
roll conditions are valid. An equation that obviously can 
be solved explicitly in the classical regime. It leads to. 



1 + 2Z.2|;^HC+^|2J 



1/2 



(47) 



where is the value of x at horizon crossing. It is im- 
portant to note here that large NG will classically develop 
when v^N,, is larger then unity whereas the self-coupling 
of X at horizon crossing is driven by . There is thus a 
domain in parameter space where the nonlinear evolution 
of the field fluctuations is essentially classical. 

Using equation (|47p and assuming that x^^ is Gaus- 
sian distributed, one can compute explicit PDFs of xi- 
They are presented on Fig. 2] for specific values of the x 
parameters. One observes that for significant values of 
the coupling constant, the large excursion values of xi 
are strongly suppressed. 

It is then tempting to assume that the resulting field 
properties can derived by the application of such a local 
transform. That would indeed lead to interesting phe- 
nomenology regarding the local metric preperties. It is 
however an approximation the importance of which is dif- 
ficult to grasp. Indeed the relation (j47|l would be valid 
if the fluctuations had all the same scale and cross the 
horizon at the same time. This is not so. For instance 
when one observes fluctuation at 10 Mpc scale, its nonlin- 
ear evolution is bound to be affected by modes at much 
smaller scales, that have crossed the horizon later, but 
that nonetheless affect the nonlinear evolution of the field 
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by changing the local field values. In a field theory lan- 
guage, this is a UV effect in the radiative correction. Its 
importance cannot be neglected since one expects modes 
that are up to e^*^ smaller than the observed ones to cross 
the horizon during the supre-Hubble evolution of the ob- 
served modes. 

The detailed description of these effects is left for a 
forthcoming study. The transform (j47p however suggests 
that the main effect of such a fully nonlinear evolution of 
the field is to squeeze the rare event tails. Such an effect 
would be visible from the behavior of the moments of the 
distributions shown on Fig. [J] 



IV. CONCLUSIONS 

We show here that it is possible to build models of 
multiple field infiation in a supersymmetric context that 
can generate a significant amount of non-Gaussianity. 
The models we advocate here are natural extension of 
D-term inflationary models. As models producing pri- 
mordial non-Gaussianities, these types of models fall in 
the same class as those described in [l^. The only dif- 
ference is that the initial isocurvature fluctuations are 
produced by a complex scalar fleld (not a real) which 
(mildly) affects its high-order correlation properties. 

Such models are characterized by a mixing angle and 
an intrinsic self-coupling parameter. In SUSY context, 
when one wants to avoid fine tuned terms in the su- 
perpotential, the potential is naturally quartic and such 
models are remarquably constrained. In particular the 
shape of the inflationary potential has only a limited 
number of parameters. We note that in such a setting 
the mass of the isocurvature modes is automatically pro- 
tected against radiative corrections (as an application 
of the results presented in [13] )■ In our analysis we 
have focused on a regime suggested by the observational 
constraints regarding cosmic strings formations. Those 
constraints impose that the Hubble scale is essentially 
constant during the inflationary phase which makes the 
transfer of modes, from isocurvature to adiabatic, poten- 
tially very efficient. Finite volume effects also introduce 
new dynamically determined parameters that induce ef- 
fective cubic terms in the potential jsij the consequence 
of which we also present. 

The phenomenological consequences of such a family 
of models can be obtained analytically. The induced high 
order correlation functions are due to super-Hubble evo- 
lution of the field. This opens the possibility of having 
weakly non-Gaussian fields. Two parameter domains can 
be distinguished for its phenomenology. If the mixing 
angle is generic then the current constraints suggest that 
the coupling parameter is small enough so that a tree or- 
der calculation of a perturbation theory suffice to derive 
the bi- and tri-spectra of the metric field. On the other 
hand if the mixing is small enough (0.1 for instance), 
then the resulting metric statistical properties might have 
been shaped by a nonlinear evolution of the field during 



its super-Hubble evolution. In principle, such an evolu- 
tion can nonetheless be addressed analytically because, 
at such scales, the fields enter a purely classical behavior. 
We present some tentative description of the evolution of 
the high order moments of the field. 

The set of predictions we have obtained has eventu- 
ally been confronted to the current constraints and to 
the expected constraints that the Planck satellite is ex- 
pected to provide us for. It is shown that a joint use 
of CMB temperature bispectrum and trispectrum is re- 
quired to efficiently explore the parameter space of such 
models. This, with the curvaton mechanism, is one of 
the very few explicit models in which primordial NG can 
significantly exceeds those naturally induced by the grav- 
itational dynamics. In this case, it is fully determined by 
the inflationary sector of the theory. 

V. CORRELATORS OF A TEST FIELD IN DE 
SITTER SPACE 

We review here the formal expression of the connected 
bi- and tri-spectrum of the component of a test scalar 
complex field in a de Sitter background. We assume that 
the potential contains cubic and quartic terms, 

V{x) = Qxt + Xlxi + XIX2X2 + ■ • • ) (48) 

as implied by the form (j32p . The dots represent terms 
that won't affect the (tree order) expressions of the bl- 
and tri-spectrum of xi- 

The perturbative calculations of those correlators can 
be done with the In-In formalism (see 2, 28]). If Q is 
the quantity the vacuum expectation value one wants 
to compute then it can be shown that a perturbative 
expansion is given by 

(o]Q(x(^))]o> - r Am r Am... 

Jrio Jrio 

m ■ ■ [[Q{X^°\^)), ~iH^'Hvi)],~-^H^'Hv2)] ■ . .]]0), (49) 

where 7f (rj) is the interaction part of the Hamiltonien, 
rjo is an arbitrarily early time when the interaction is 
supposed to start playing a role, r] is the time at which 
the expectation value is computed and x^^"^ iv) ^re the 
field values at time ry when they evolve according to 
the non-interacting part of the Lagrangian. In the fol- 
lowing we use this formulae for Q = xi(ki) . . .Xi(k3), 
Q = xi(ki)---Xi(k4) and for ifW = Jd^^y^Vix) 
while keeping the expansion only to tree order. 

The expression of those correlators will depend on the 
time dependance of the free field. The latter is given by 

xf )(x) - ^ / d^k (A(77)ake"'- + mv)ale-''^-) (50) 

where, 



10 



for a de Sitter background, rj is here the conformal time, 
77 = — 1/ (aH) . this function permits to define the differ- 
ent time Green function, Gfe(r/,7/), 



Gk{r],v') = fk{v')fk{v)- 



(52) 



After some calculations, one can show that the tree 
order term for the three-point function is given by. 



r" dr]' 



1 



J [Gk, iv, v')Gk, iv, v')Gk, (ry, ry') - c.c] (53) 



-00 Hr]' 



It is interesting to define a reduced bispectrum 
Q3(ki,k2,k3) as 



Q3(kl,k2,k3) 



P3^^(ki,k2,k3) 



Pxi(fci)pxi(fc2) -t-sym. 



(54) 



This function is shown on Fig. [5] as a function of ry, 
or more precisely as a function of the number of efolds 
since (or before) horizon crossing. The mode configu- 
ration correspond to an equilateral configuration (ki = 



^2 = ^3) but the asymptotic behavior of Q3 is indepen- 
dent on the configuration (dashed line) if is defined 
as A^e = — log(— fctjy) with kt = ki + k2 + ks. 

For the four point function two types of diagrams are 
contributing at tree order. One is due to the quartic 
part of the potential. In the following we will denote it 
as the " star" contribution. The other is due to the cubic 
terms of the potential. There are two such contributions 
(one due to an exchange of xi degree of freedom and one 
to an exchange of X2)- Altough a priori smaller than the 
previous term, these two contributions are not necessarily 
negligeable. 

The star contribution is given by [3^ . 



Pg'^4j.(ki,k2,k3,k4) = -6iz.2 



dr^' [Gfe, (77, ry')Gfe, (ry, r^')Gk, (77, v')Gk, (ry, ry') - c{6^) 



The line contribution is finally given by the formal ex- 
pression. 



J 



Pj^;^g(ki,k2,k3,k4) - -2,y^{9xl+xl) 



dryi d?y2 
Hrji Hr]2 



[G|ki+k2|('yi,'?2) +C.C.] [Gk,{v,Vi)Gk2{v,m) - c.c] [Gfc3(7y,7yi)Gfe4 (77,772) - c.c] -1-sym.. (56) 



Similarly to Q3 one can define the reduced tri-spectra This function is shown on right panel of Fig. [5] for a pecu- 
as[42|. liar configuration. Note that the relative importance of 

the two depends on the configuration even in the asymp- 
totic limits. 



'3star(ki,k2,k3,k4 



n?ar(ki>k2,k3,k4 



Pxi (fci)Pxi (yc2)Pxi (fcg) + sym. 



(57) 



and 



'3line('^i'^2,k3,k4 



P-^(k„k2,k3,k4 



Pxi(fci)Pxi(fc3)Pxi(|ki +k2)| +sym. 



The asymptotic limit of the line part of the tri- 
(58) spectrum is given by, 

I 



Pji^e (ki , k2 , k3 , k4 ) ^ - _^ „4 , 3 13 Jla i.3 [94 (ki , k2 , 77)54 (k3 , k4 , 77) sym.] (59) 

with 

(74(ki,k4,7y) = 2(fci+fc2)fc2 + (fc3 + fc3)(2-27) 

+2kik2{kl + kl) fc3log[(fci + fc2 + k)/{ki + k2- k)] - {kl + ki)log[r,^{ki + k^ + k){ki + k2 - k% (60) 

I 

where /c — |ki-|-k2| and 7 is the Euler number. Note that the dominant contribution of at super-horizon scale 
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equilaral triangle configuration square configuration 




-2. 0. 2. 4. -2. 0. 2. 4. 

-Log[-kt I/] -Log[-kt 17] 



FIG. 5: Behavior of the reduced correlateors, Q3 and Q4 of the field xi; as a function of iVe = \og{ktri). The function Q3 

(left panel) is to be multiplied by u'^Xi/H'^^ the function Qf^^'^ (thick lines of right panel) is to be multiplied by —v^/H^ and 

the function Q^i^^ (thin lines) is to be multiplied by 2i/''(9xi +X2)/-'^*- The dashed lines correspond to the corresponding 
asymptotic behahiors. 



is due to the last term. It behaves like —2{kf + fc|)A^e- 
This leads to the super-horizon behavior of the four point 
function of xi exploited in the text. 
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